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ABSTRACT 


The validity of linear smoothing of topography for numerical weather 
prediction and the variation of mountain drag with mountain height and 
Static stability are examined in this study. In the model a constant 
geostrophic current is perpendicular to the mountain range and the 
height of the mountain is independent of y. The hydrostatic Boussinesgq 
equations are used with motion bounded at the top by a rigid plane at 
z =D. A modified coordinate system similar to Phillips’ sigma system 
was used. Solutions were obtained using a smoothed mountain profile. 
These solutions were compared with smoothed solutions obtained from the 
unsmoothed mountain. The comparison of these solutions shows that an 
error is introduced when non-linear terms become sufficiently large. 
Values of the mountain drag for differing values of mountain height at 
a given static stability and for differing va lues of static stability 
at a given mountain height were computed. Mountain drag was found to 
vary quadratically with mountain height and linearly with static 


stability. 
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I. INTRODUCTION 


The use of large scale numerical forecast models makes it necessary 
to introduce a topographical model of the surface of the earth for use 
in the lower boundary condition. Cressman [1960] notes that since the 
topography of the surface of the earth has many small scale features, 
it is necessary to smooth the topography when it is used as a boundary 
condition in a numerical forecast model. Mintz [1956] and Cressman 
[1960] indicate that the influence of small scale features of the topog- 
raphy can be included in a surface drag coefficient which varies from 
point to point. 

One purpose of this study is to examine the validity of linear 
smoothing of the topography in numerical weather prediction. To accom- 
plish this the flow over a small scale mountain was computed using a 
small grid mesh. The results of these computations were then smoothed. 
The mountain was then smoothed and the solution was recomputed. Since 
the parameters predicted in large-scale forecasts are interpreted as 
area averages, a comparison of the two results should show little, if 
any, difference if the smoothing technique is valid. 

Mountain drag, or mountain torque, which arises from a difference 
in pressure on the two sides of a mountain, is important in the general 
circulation of the atmosphere. Studies of mountain drag associated with 
smaller scale mountains could be used to improve the calculation of the 
surface drag coefficient which was mentioned above. Miles [1968] showed 
that wave drag on an obstacle is a function of the height of the obstacle 
and of the mean flow perpendicular to the obstacle. The present paper 
investigates the effect of mountain height and initial static stability 


on mountain drag. 
Jd. 








This study employes a two-dimensional, primitive equation model. 
The geostrophic current perpendicular to the mountain is constant. 
The mountain range is infinite in lateral extent which makes the model 
independent of y. A special coordinate system similar to that developed 
by Phillips [1957] was used. The equations were solved numerically 
with the use of a special finite difference scheme which was deve loped 


by A. Arakawa (see Langlois and Kwok, [1968]). 
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II. THE BASIC EQUATIONS 


‘The basic equations used in this study are the Boussinesq equations. 
These equations will be applied to a model in which the flow is bounded 
above by a rigid horizontal plane. These equations neglect the compres- 
sibility of the atmosphere which, for this study, should not be of 
qualitative importance. It will be seen that the internal scale of the 
flow is small compared with the depth of the fluid. The hydrostatic 
approximation is made and heating and friction terms are neglected. 

The Boussinesq equations (Ogura and Phillips [1962]) may be expressed 


in the z-coordinate system as 


ee eee (1) 

2+ Fverw =o, (2) 

Sk og, (3) 
: 

vi + 2 = 0 (4) 


where: 


9g = T(p,/p) ~ Q. 


In this model the only quantity which varies with latitude is $ and it 
varies linearly. The Coriolis parameter is constant. 

To simplify the finite differencing of the basic equations a 
modified coordinate system is introduced in which the lower boundary 
surface becomes a coordinate. This system makes it possible to intro- 
duce the effects of topography without employing uncentered horizontal 


space differences in the vicinity of the mountains. [It also allows 


Ls 





application of the lower boundary condition at the actual surface of 

the mountain. In the proposed coordinate system the vertical advection 
terms are identically zero at both the upper and lower boundaries. This 
system of coordinates is very similar to the sigma system introduced by 
Pivilips (1957 |. 


The vertical coordinate, zeta, is defined as: 


Z- z 
5 (5) 
where: 
Z=D- 2 
S 


According to this definition €¢ = 0 at the lower boundary and € = 1 at 
the upper boundary. Fig. 1 illustrates this coordinate system. 
Transformation of equations (1) through (4) from the z-coordinate 


system into the C-coordinate system yields: 


Tega we % (C= 1 2 i 
-ve- F4Es*)ye,- kx (6) 


2+ tve+¢ V=0 (7) 
a oe (8) 
a6 8. 

vat 2-0 (9) 


In equations (6) through (9) V is the horizontal gradient operator on 
zeta surfaces. 

Since some of the following computations involve vertical averaging 
it is convenient to express the equations in flux form. Equations (6) 


through (9) in component form become: 


OM) 
x tes - 28th S+ ay (10) 
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MH, ) 

St sus uae wee as 
OH | _ 

St + © (8) = 0 (12) 
of _ gs 
x 7 6 H (13) 
OM) ig 
Saas + Z xX = 0 (14) 


where: 


e(s) =< (su) + $= (¢s2) 


and M = uZ, M_ = vZ, H = @Z. 
x x 


The boundary conditions are: 


C =O at C¢ 


C 


O, and (15) 


One ee1G Le (16) 
For this study the flow will be independent of y and periodic in x. 
Since ¢ is the only variable allowed to change in the y direction 

x = 0 for all other variables. 


To obtain an expression for $ equation (13) may be integrated with 


respect to C, 


O 


C 
e-(2 nate (17) 
0 


The constant of integration, C, is eliminated by taking the vertical 
average of equation (17) and subtracting which gives the following 


expression for #4: 
c 


: 
eh eee ae) * (18) 
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In the above equation and throughout the remainder of this paper | 


represents the vertical average of a given variable and — represents 
the, horizontal average. 

| 

' Before a solution for % can be found, a value for ®@ must be obtained. 
This is accomplished by first taking the vertical average of equation (10) 


thereby eliminating the vertical advection term. Next the vertically 


averaged equation is differentiated with respect to x yielding: 


2 2 a 
o rom, 9d My el OZ rv 
Se (M,) + wae 2 |- la - x (C-1) = or va (19) 


The first term of this equation is zero since the vertical average of 


equation (14) vanishes and the resulting equation is: 


.¢O-bik a+ Pecos ml 
a Z Ee (M ) + H (C-1) Ss gee (20) 
This equation can be solved diagnostically for 9. 


To obtain a equation (18) is differentiated with respect to y 


ylelding: 


Od ; 00 ; ag ory 
= [ehia e Oe 
= o Z Ys dc = ac) + Sy 21) 
Since <2 = 0 for this study 
0d _ J? i (22) 


and this quantity is assumed to be independent of x to keep the 
perturbations independent of y. 
Now consider the domain averages of equations (10) and (11) which 


abe; 


d_ _ Bo OZ 
dt e728 ey (23) 
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@) = - 2 = - ‘mM (24) 


CM) (25) 


If this relation is used for all time and if the mass flux in the y 
direction is zero initially, then © will be zero for all time. The 
first two terms on the right side of equation (23) represent the mountain 
drag. The sum of these two terms was shown by computation to be negative 
in the mean as expected from physical considerations. This will lead to 
a continual loss of x-momentum. A loss such as this would make a quasi- 
steady final state impossible. An artificial source of momentum is thus 
required in this simple model to maintain the mean zonal flow. The 


1 


artificial source of momentum is provided by holding a as obtained in 


equation (25), constant for all time. In this case if Mis decreased 
below the initial geostrophic value MY will increase and restore the 


momentum loss. Thus, in the long term mean, the mountain terms are 


balanced by the mean Coriolis force due to a net y-momentum flux. 
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IIil. FINITE DIFFERENCE EQUATIONS 


Equations (10) through (14), and (18) are solved by introducing 
finite differences in x, €, and t. The space finite differences employ 
a variation of a scheme developed by Arakawa [1966] (see also Lilly 
{1965]). The basic scheme was designed to conserve total energy in 
addition to conserving the squares of u, v, and 6 with respect to the 
advection terms. This scheme was also designed to avoid the non-linear 
computational instability which was discovered by Phillips [1959]. 

The time differencing scheme employs a variation of the Matsuno 
[1966] two stage differencing scheme which is also referred to as the 
Euler simulated backward time step differencing scheme. This scheme is 
similar to that used in the general circulation model developed by A. 
Arakawa (see Langlois and Kwok [1968]). The modification makes it 
possible to avoid a checkerboard spatial pattern by alternating between 
centered and uncentered space differences in certain terms. 

The Euler backward difference scheme may be written in the following 
form 


s*¥ =s + AtF® (step 1) 


S = 5° + ATF* (step 2) 
where 
Ss =u, v, 96 and 
F is the right hand side minus the linear 
operator £(u) as defined by equations (10) 
Ehisoueh (12). 
Centered space differences are employed to compute Fe in step 1 above 


and uncentered space differences, using the values of s* obtained in 


step 1, are used to compute F* in step 2. During computation, forward 
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spatial differences are used for even time steps and backward spatial 
differences are used for odd time steps. This method allows geostrophic 
adjustment to occur at the smallest scale recognized by the model. 

‘The vertical domain is divided into 18 layers such that AC = 1/18 
for each layer. The boundary surfaces are denoted by the subscript i, 
where i = 1 at the lower boundary, the surface, to i = 19 at the top of 
the model. The quantities M, Me H, and $ are computed on the even 
levels and C is computed on the odd levels. Throughout the remainder 
of this study j denotes the horizontal grid index with Ax as the grid 
mesh length. For the remainder of this paper CDF, FDF and BDF will 
represent centered difference form, forward difference form and back- 
ward difference form, respectively. 

The finite difference form of the linear operator (s), where is is 
any dependent variable, can be written | 


a 7 Kee 2a Mi) - Sap oe eel an, >| 
M(s)i 7 2 Ax " 


2; (Fe "itl, j i so (26) 
2AC 


In equation (26) the flux forms Ms and M are defined as: 


> ¢ eis ea. 
xX-+ 2 e 
(M . .+M ..,) 
= xe] x i,j-l 
Me 9 (CDF ) 
hs = ae i,jt+l ote 
M =M_. ., (FDF) and 
X- >a Ue 
X+ ~ Lt i) ist 
Mo = Mg 5.1 (BDF) 
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The finite difference approximation to the integral in equation 


Giz ) Sis: 


H dC = ; Gy) a) oe (27) 


The value of the constant of integration, C, is not important since it 
will vanish when the vertical average of equation (27) is taken. 


When i is odd, equation (13) can be approximated with: 


= g_ 


The continuity equation, equation (14), can be approximated with: 


e . e uy 2G a ; ; 
Gi41,j a Sets as Z Ax (Ma M_) (29) 


The flux forms Mo 


Since all of the dependent variables occur at each grid point, the 


and M_ are as defined above. 


x derivatives are approximated by: 


®, 8 a0 


1,J 


The form of the dependent variable §s is defined as: 


Crete 
ee ee (CDF) 


6s = enn = Sip (FDF) and 


Ss (s. .°7 S. fet (BDF ) 


1,J 


To maintain consistency within the model the finite difference form 
of equation (20) is obtained from the finite difference form of equation 


(12) and appears in the following form: 
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a - . oD 
A em Re (31) 





where: 
B. =A, +C 
J d J 
A, = Z. and C. = Z. used with CDF 
j jtl j ee) ) 
A. = Z. and C, = Z, (used with FDF) 
j j j ao 
A. =Z. and C. = Z. (used with BDF) and 
j aes i‘ 
PE a a rey | 
p, = A (-iettl had - id Siat-l 
j Ax 2 x Z x- 
g 5-7 6Z I 
+ ., = + MM .. 
eis Ax y i,j 
where: 


(D.., - Dy_,) 
6D = tt (used with CDF) 


6D =D. - el (used with FDF), and 


6D 


Il 
oO 


ee D. (used with BDF) 
In the above equations M.. , M_> and §s are as defined above. 


<< 
A solution for $ can now be obtained using the exact technique 


described on page 198 of Richtmyer and Morton [1967]. 
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IV. SCALE ANALYSIS 


The following scale analysis is done in the z-coordinate system as 
a matter of convenience. The analysis will be useful in determining 
the initial conditions for this model and will serve as an aid in the 
interpretation of the numerical solutions. 


The following non-dimensional independent variables are defined: 


x' = a and (32) 
C Se (33) 


L and d are the horizontal and vertical scales of the disturbance. 


They are due to the characteristics of the mountain and other physical 


parameters. 
u=U+R* vu! (34) 
v= A. Lfv' (35) 
t 
we i) (36) 
30. 
on ee | (37) 
a a 30. 
= — Zz! 
$ f°, dz + Bd Q- x 6 (38) 


ple constant, U, Tepresents the zonal flow away from the mountain. The 
5 S cee é ; : 

quantity oz > represents the stratification of the basic undisturbed 

atmosphere which, for this study, is assumed to be constant. The scale 


of the meridional component of the wind is: 
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Two dimensionless numbers appear in the analysis. They are; 


m UD 
R, = FL and (40) 
Fy Yee 
Ry ~ £1, Seq (41) 


Ro is the Rossby number based on the zonal wind and the mountain scale. 
Ro is the disturbance Rossby number and is so Called because it is based 
on the amplitude of the v component of the wind. 


The scale analysis requires that: 
a 30. ue 
d = Lf (2 Seu ) (42) 
O 
as long as d < p. 


Using the above definitions the steady state non-dimensional 


equations may be written: 


Ry a ur Se wi SE) a Be ys (43) 
srt 8, (w oa ch Sy) tu! (44) 
Soo wh ER ut SO- 4 ys 20" = ( (45) 
sr + Se : (46) 
sr a (47) 


dz! 
w! 2 ea at z = z, and (48) 
w' = 0 at z! =? (49) 


24 





25 


> terms are ¢liminateg , the 


If just the R* Cerm is dropped , 


Tney [1962] result, 





V. INITIAL CONDITIONS 


The profile of the mountain used in this study is given by the 
cosine squared of the horizontal distance. The mountain is 1200 km 
wide, has an amplitude designated B and is duplicated every 6000 km. 
The lower boundary of the model is flat and horizontal except where the 
mountain exists. 

The initial state is derived from the steady state, quasi-geostrophic 
equations which were shown to be linear in the preceding section. These 


equations are: 


yee 
al Ox ! (50) 
u’ = - Ovi (51) 
So" + wt = 0 (52) 
Our. Owe.) 
ox ae (53) 
ob’, 
S50 wa (54) 


The equation necessary to obtain the initial fields of u, v, and 
8@ is the potential vorticity equation for this set of linearized 


equations. To obtain the potential vorticity equation, differentiate 


4 y ° 
Di a 


equations om and (52) with respect—to-x. Subtract equation (53) 
from the differentiated form of equation (52). The resulting equation 


is: 


aE) Be 09 
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Now differentiate equations (51) and (52) with respect to x and equation 
(54) with respect to z and substitute into equation (55). The resulting 


equation is the potential vorticity equation: 





2 2 
0 ono Ce 
set Sa + Sa)" ; ti 


To obtain the initial solution for this study the following boundary 


conditions are applied: 


| 

= =Oatz' = > and (57) 
| 

5 =z, at z= 0 (58) 

Equation (56) is satisfied if: 
Jit ee 
oz! ox' 
A general solution of equation (58), here x' = kx, is: 
6" (x'\2') = F(X') cos x! (59) 


Applying this solution to equation (58) at the upper boundary and 
solving for F(z') yields: 
iNet ' D 

F(z‘) = A cosh (2 - 2’ (60) 

From this it follows that: 
S | | | | D 
$6’ (x',z') =A cos x cosh (2 - 2) (61) 
The mountain topography can be Fourier analyzed into a cosine 

series. Each component of that series will have the following form 


if x is non-dimensionalized with the wave number of the component. 


Zo = B cos x (62) 
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When the solution for this component is obtained, the quantities can 


be dimensionalized and summed to satisfy the complete boundary condition 
‘Differentiation of equation (61) with respect to z' and evaluation 
' gives the following result: 


of the resulting equation at z' =z - 


i] 
SU) = A cos x' sinh ('. = >) (63) 
S 


If it is assumed for this analysis that 2 is sufficiently small that 
0, then this equation can be 


equation (63) can be applied at Zz" 


solved for A: 


A=- ee (64) 
: D 
sinh (2) 
The final solution for ¢' is: 
(65) 


G' = cos x' cosh GS - >’) 
sinh (3) 


From this solution for $' follow the solutions for u', v', and Q': 


vi=- — sin x' cosh (2 - 2) (66) 
sinh (- ) 
i = ee cos x' cosh (2 - 2) (67) 
: D 
sinh (3) 
(68) 


Q' = —— cos x’ sinh (2 - 2) 
sinh (2) 
A Fourier series of the full lower surface yields: 


z = 2% B(k) cos kx 
= nk 


in dimensional form. Solutions (66) through (68), when dimensionalized, 


can be summed to give the following initial conditions: 
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Q@ dz 
i= a = ee [a(z-D)] cos (kx) (69) 
k=1 £° q sinh (aD) 
a Pee 
nw Be TELS) oe Oz 
v= aa f sinh (aD) cosh fs sin (kx) (70) 
ae 
nw B(k) — 
@=@0 + 2& Seen) sinh [a(z-D)] cos (kx), and (71) 


B sin (™ npm T) 
B(k) = eee (72) 


ow [1 - 22) ] 


where: 


nw is the number of wave numbers used, npm is the number of points in 
the mountain and np is the total number of horizontal grid points in 
the model. 

The first 20 wave numbers were used to initialize the wind and 
temperature fields. Any number up to and including 20 could have been 
used with k = 20 yielding the smoothest initial solution. Any number 
greater than 20 would have introduced aliasing since the waves generated 
would have a wave length less than two grid lengths. 


Me -l 
In the initial fields the constants used were U = 8m sec , 


08. O -1 -2 ,o,,-1 -4 
sy = 4.0 “K km, D = 9 km, = = 0.0327 msec (°K), and f = 10 

O 
sec . The value of 0. was computed for each point from the relation: 


32 [Pe2 G0) 
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VI. RESULTS 


In this study Ax = 150 km, At = 20 min, and pe equally spaced 
levels in zeta ane used. Seven values of B vith = = 4 °K elie and 
eleven values 34 at B = 1575 m were employed. 

Values of the mountain height were computed using the function 
B eos” (kx). This function provides a smooth profile at the base of 
the mountain. First, solutions none computed using this function to 
determine the mountain profile. The solutions were then computed using 
a linearly smoothed mountain profile. The smoothing was accomplished 


by applying the following linear formula: 


Vv 
+ 2. + 2 
VEL eu Re > 


This procedure was repeated four times. The resulting mountain height 
was about two thirds that of the original mountain height and the 
horizontal extent was approximately doubled. In each case the solutions 
were allowed to converge toward a quasi-steady state. To permit com- 
parison of the various cases the results were converted from the zeta- 
surfaces to the z-surfaces at one km intervals. In this study R is 
independent of the mountain height and its value is RS = .419., 

The general features of the flow were basically the same for all 
cases and were similar to those reported by Walton [1968]. There was 
a ridge in the pressure field over the mountain. The v component of the 
wind was from the south on the windward side of the mountain and from 
the north on the leeward side. The u component of the wind was at a 
maximum near the surface over the mountain and decreased rapidly with 


height to a minimum at the top of the atmosphere. Over the plain, 
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away from the mountain, the u component was a minimum near the surface 
and increased very slightly with height to a maximum at the top of the 
BEmoS pucker The potential temperature was a minimum over the mountain 
aoe 

Fig. 2 is a graph of the difference between the maximum value of $ 
in the ridge and the minimum value of ¢ at the bottom of the trough, at 
the three km level, with respect to mountain height. The upper curve 
represents the smoothed solutions from the unsmoothed mountain. 

The difference between the smoothed and unsmoothed solutions is 
least for mountain heights less than 900 m. Below this height the 
solutions are almost identical. The scale analysis of the basic 
equations used in this study showed that when R* is small the balance 
approximation should hold. When R is also small the equations reduce 
to the linear quasi-geostrophic equations. Since the equations are 
linear for ’R | small the solutions should be nearly equal in this range 
as they actually are. At mountain heights greater than 900 m the 
non-linear effects become large enough to cause a marked separation 
of the two curves. These results are comparable in magnitude to those 
obtained by Walton [1968]. In his study he applied a different differ- 
encing scheme to obtain solutions to this same model. 

The addition of the artificial source of momentum suggested in Section 
II appears to have been valid. The slow decay in mean zonal velocity 
with time noted by Walton [1968] did not occur. In providing the 
artificial source of momentum it was assumed that the model would come 
into adjustment and that the mountain drag term would be balanced by 
the term denoting mass flux in the y-direction. Fig. 3 shows a plot of 


these two quantities versus time. This figure shows the values for a 


a1 
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single mountain height, however it is representative of all the mountain 
heights in this study. The upper curve represents the momentum flux in 
the | y-direction and the lower curve represents the mountain drag. It 
Barcus seen that these two quantities do behave as expected after the 
initial period when the model is coming into adjustment. It should be 
noted that the sum of the two curves does not equal zero. The long 

term change in the x-mass flux is very small. The reason for this 
difference is apparently due to the fact that centered space differences 
were used to calculate the mountain drag while uncentered differences 
were used to calculate the actual forecast variables. 

The solutions, represented by A?, did not appear to be converging 
significantly at the end of 40 days into the forecast period. For the 
higher values of mountain height there appeared to be slight conver- 
gence, however this may have been due to the time finite differencing 
scheme which tends to damp the computational mode of solution. Fig. 4 
is an example of the oscillation of the solution obtained. These 
oscillations have a period of about eight days. This pattern was 
identical for all cases with only the amplitude of oscillation varying. 
The amplitude increased with increasing mountain height. This oscillation 
appears to be the result of mechanical interaction between the mountains. 
One case involving a mountain height of 2025 m was computed retaining 
the width of the mountain at 1200 km but doubling the width of the plain 
from 4800 km to 9600 km. As a result of this the number of oscillations 
occurring in a 40 day forecast was diminished by one half with the ampli- 
tude of oscillation decreasing slightly. The values used in Fig. 2 are 
the mean values, at a given height, about which Aé¢ oscillated. 

As a check of the validity of the geostrophic assumption the v com- 


ponent of the geostrophic wind was computed. Except in the near vicinity 
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of the mountain the values obtained from the forecast equations and the 
values obtained from the geostrophic computations were similar. The 
difference between the geostrophic and actual v component was less than 
ten percent in conformity with equation (13). 

The variation of mountain drag with mountain height and with 
static stability was also investigated. Fig. 5 displays the results 
of the variation of mountain drag with mountain height. From this 
figure it can be seen that mountain drag increases quadratically with 
increasing mountain height. Fig. 6 shows the relationship between 
mountain drag and static stability. This figure shows that mountain 
drag increases linearly with increasing static stability. An investi- 
gation of the gross features of the flow for increasing static stability 
at a constant mountain height showed no significant changes except in 
magnitude. 

The stream function was computed and streamlines were drawn for 
each mountain height at selected time spacing. Fig. 7 is a typical 
example of the streamlines obtained. It shows horizontal flow over 
the plain at all levels and flow conforming to the mountain in the 
vicinity of the mountain with the upper level flow nearly horizontal 
everywhere. 

In the case where B = 2250 m static instability developed about 
36 hours into the forecast period. A comparison of the u and ¢$ fields 
obtained in the hydraulic jump case B described by Houghton and 
Kasahara [1968] with those obtained in this study for B = 2250 m shows 
a marked similarity in the structure of the fields. Houghton and 
Kasahara used a single layer, non-rotating model. The lower values of 


mountain height used in the present study apparently fall into the 
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non-critical, or no jump domain, described by Houghton and Kasahara. 
As mountain height increases the values eventually fall in the hydraulic 
jump domain. The reason for the static instability at B = 2250 m is 
Bits attributed to hydraulic jump. Fig. 8 shows the values of u 
initially and six hours later for B = 2250 m at the three km level. 
The development of the minimum in u in the lee of the mountain is 
already evident even though static instability has not yet occurred. 
This model is not appropriate for the treatment of hydraulic jump 
because of approximations which have been made. Calculations of 
mountain drag were not made for this case but it is expected that if 
these calculations had been made the character of the drag properties 


would have been modified above B = 2025 m. 
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VII. CONC LUSIONS 


Since this study involved only one horizontal scale, only one 
yalue of the Rossby number (R, = U/£L) was considered. Solutions 
obtained in this study showed that linear smoothing of the topography 


an error which increases as the disturbance Rossby Number 
14 


[z, == aus ° —=) | increases. If the Seratiies wee aie study 
could be extrapolated to smaller scales, i-e-> where Ke is larger, the 
difference between results obtained using 4 smoothed topography and 
those obtained using an unsmoothed topography would be even greater. 
This is consistent with the example given by Charney (1967] in which 

he considered an infinitesimally thin mountain, oriented along 4 
meridian and extending to the top of the atmosphere. Without smoothing 
it would act as an absolute barrier, however, with smoothing its 


vertical extent would be decreased to nothing and it would have no 


effect on atmospheric flow. However it is not clear what error will 


—_ 
i] 


be found if RE becomes large enough for pressure jumps to form. The 
error behavior will probably be different than that obtained here. 

The addition of an artificial source of momentum did maintain the 
mean value of the u component of ae flow at a value near the initial 
mean zonal flow. A comparison of the values of R, obtained in this 
study and those obtained by Walton [1968] did not show any significant 
difference in the results obtained for smooth and unsmoothed mountains. 

By Aouigline the size of the plain between the mountains it was 
observed that the period of esc ilatrons on cue Ad field was reduced by 
one half. It was thus concluded that this oscillation is due, at least 


in part, to mechanical interactions between the mountains. 


42 : 





Mountain drag computed for varying mountain heights at a constant 
Static stability showed an apparent quadratic relationship between 
increasing drag and increasing mountain height. Mountain drag computed 
for varying static stability at a constant mountain height showed an 
apparent linear relationship between increasing drag and increasing 
stability. 

For the case when B = 2250 m static instability resulted early in 
the forecast period. The reason for this instability was attributed 
to the formation of a hydraulic jump. This cannot be definitely stated 
as the model, in the form which was used for this study, was not 
appropriate for the treatment of hydraulic jump because of the lack of 
friction and because of certain approximations which had been made. 

Further studies with this model should be conducted using a wider 
range of mountain scales. This would give a more complete relation- 
ship between mountain drag and mountain structure which could be used 
in parameterizing mountain effects in numerical weather prediction 
models. The drag could be obtained for the larger values of KR. if 
appropriate treatment of the hydraulic jump is included. The mechanical 
interaction between the periodic mountains could be reduced if a variable 
mesh length grid in the horizontal scale of the model is used. This 
would permit a short mesh length in the vicinity of the mountain where 


most of the activity is. 
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